Melting a granular glass by cooling 
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Driven granular systems readily form glassy phases at high particle volume fractions and low 
driving amplitudes. We use computer simulations of a driven granular glass to evidence a re-entrance 
melting transition into a fluid state, which, contrary to intuition, occurs by reducing the amplitude 
of the driving. This transition is accompanied by anomalous particle dynamics and super-diffusive 
behavior on intermediate time-scales. We highlight the special role played by frictional interactions, 
which help particles to escape their glassy cages. Such an effect is in striking contrast to what 
friction is expected to do: reduce particle mobility by making them stick. 

PACS numbers: 



Friction implies stability. A solid block only remains 
static on an inclined plane because there is friction with 
the underlying surface. Similarly, a heap made of granu- 
lar particles, in general, will have a higher angle of stabil- 
ity with increased friction coefficient [l| . Inside the heap, 
the particle volume fraction will generally be lower, and 
the number of inter-particle contacts smaller, while the 
heap nevertheless remains stable 0. A similar stabiliza- 
tion occurs when the grains are driven into a fluid state. 
Under shear, the jamming transition from a freely flow- 
ing state to a yield-stress fluid occurs at lower volume 
fractions as compared to the frictionless case 

Here, we present simulations of a driven granular sys- 
tem were friction acts opposite to what is expected from 
these simple examples. We show how friction can lead 
to anomalous particle dynamics that very efficiently flu- 
idize the system. As a result the system undergoes a 
re-entrance melting transition from a glassy to a fluid 
state by lowering the amplitude of driving. 

Experimentally, a variety of driving mechanisms have 
been proposed to characterize the dynamical properties 
of dense granular systems. Among those are shaking Q 
fluid- or air-flow H -fuji cyclic shear [ll], [l2| or tempera- 
ture oscillations jl3l . Il4j . We describe a two-dimensional 
system, similar to the setup used in Refs. [H, 0, [13, EE)]. 
At high densities and low enough driving amplitude these 
systems readily form glassy states, where structural re- 
laxation is completely suppressed flij . Interestingly, 
Ref. [HI also reports anomalous particle dynamics deep 
in the glassy phase, and suspects friction to play a central 
role in this process. Unfortunately, it is rather difficult 
experimentally to quantitatively characterize or tune the 
frictional interactions between the particles [l7j] or be- 
tween particles and container walls. The connection be- 
tween friction and particle dynamics is therefore unclear. 
Our simulations have the goal to elucidate such a con- 
nection. 

Model - We consider a two-dimensional system of 
N = 2500 disks. One half of the particles ("small") 
have radius R s = 0.5d and mass m s = p(4n/3)B%, the 
other half ("large") have radius Ri = 0.7d and mass 



mi = p(47r/3)i?f. Units are chosen such that particle 
mass density p = 1 as well as diameter d = 1. The 
particle volume fraction is defined as <fi — Yli=i ^Ri/L 2 , 
where L is the size of the simulation box. Unless other- 
wise stated we fix <j) = 0.825, which is only slightly below 
the random close packing value of 4> c = 0.84. Periodic 
boundary conditions are used in both directions. 

Particles interact via a standard spring-dashpot inter- 
action (e.g. 0, 0, i). In short, two particles i,j in- 
teract when they are in contact, i.e. when their mu- 
tual distance r is smaller than the sum of their radii 
Ri + Rj . The interaction force has both a normal com- 
ponent F n = k n (r — (Ri + Rj)) and a tangential compo- 
nent F t = k t S t , where k n = 1000 and k t = 2fc„/7 are the 
spring constants and 6t is the tangential (shear) displace- 
ment since the formation of the contact. The tangential 
spring mimics sticking of the two particles due to dry 
friction. The frictional forces are limited by the Coulomb 
condition F t < pF n , which is implemented by rescaling 
the tangential displacement 5t —y p,F n /kt whenever nec- 
essary. With these interactions we perform molecular dy- 
namics simulations using LAMMPS [3] and a time-step 
of At = 0.001. 

To mimic the presence of an external container, we 
place the particles on a flat surface (sy-plane) , which 
has similar frictional properties as the particles. The 
only distinction is that the particle-surface normal in- 
tcraction is set by gravity, F r \ — mig. As a conse- 
quence, the surface exerts a permanent frictional force 
on the particles, which is limited by the Coulomb thresh- 
old, < psiriig. Note, that we do not allow particles 
to move away from the surface (no motion in ^-direction) . 

Particles are driven with an oscillating force F(t) = 
Asm(u>t) that acts (in the plane of the surface) along the 
y-direction. We use the period of the driving as our unit 
of time and set T = 2it/oj = 1. Even though this driving 
is uni-directional, we find that due to the dense packing 
the system remains roughly isotropic. 

With this kind of driving force, an isolated parti- 
cle on a frictionless surface oscillates at an amplitude 
Do = A/mu! 2 . In a dense assembly this leads to local 
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frustration as particles with smaller masses tend to move 
faster. Rearrangements result which, at the high den- 
sities under consideration, may or may not be able to 
lead to structural relaxation. It is this glassy dynamics 
that we are interested in, with the driving amplitude A 
playing the role of thermal temperature. Note, that this 
driving mechanism injects energy directly into the bulk. 
System-size is therefore not an issue and the system is 
spatially homogeneous. 

The simplest quantity which is measured from particle 
displacements is the mean-squared displacement (MSD), 

A2 (*) = ^E^(*)-^(°)] 2 ) w 

where Xi(t) — Xi(0) is the displacement of particle i in 
the direction transverse to the driving. Snapshots of the 
system are taken after every full force cycle. Time is 
therefore restricted to t = t n = nT and n integer. With 
this definition the MSD is zero when particle motion dur- 
ing cycles is periodic. 

Results - Let us first consider the case where there is 
no inter-particle friction (p = 0) and only particle-surface 
friction (fi s = 1). In Fig. Q] wc display the evolution of 
the MSD for various driving amplitudes A. 

At relatively high amplitudes (A > 3) we observe typ- 
ical signs of glassy dynamics; an extended plateau in the 
MSD, caging of particles and hopping events that eventu- 
ally lead to particle diffusion. As expected the MSD de- 
creases with decreasing driving amplitude. Surprisingly, 
this trend does not persist when the amplitude is fur- 
ther reduced. Instead, for lower amplitudes we observe a 
greatly enhanced MSD and a total dissapearance of the 
plateau. Thus, the system is fluidized, very unintuitively, 
by reducing the amplitude of agitation (also see supple- 
mentary movies). This re-entrance melting transition is 
our main result. 

The fluidization is particularly evident in the long-time 
diffusivity D := lim^oo A 2 (t)/t (see Fig[2]), which has a 
minimum in the glassy regime |22j but then strongly in- 
creases up to a sharp maximum at A* « 1. At the same 
time, the "granular temperature" T, defined as the av- 
erage kinetic energy in the direction transverse to the 
drive, is completely normal and monotonously decreases 
with decreasing A. A similar trend is visible in the MSD 
at short times. Anomalous dynamics only develops at in- 
termediate time-scales. Associated with the maximum of 
D(A*), we observe an intermediate super-diffusive regime 
that spans one to two orders of magnitude in time. 

Discussion - Clearly, frictional interactions with the 
surface are important for the observed anomalous dif- 
fusion. If we switch off these forces, then we only see 
the cage- like dynamics characteristic of the glassy regime. 
Frictional forces, on the other hand, are expected to be- 
come relevant when the driving force competes with the 
Coulomb threshold, A ^. s mg. Intuitively, one would 
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FIG. 1: Mean-square displacement A 2 (t) for different driving 
amplitudes A (arrows indicate decreasing amplitude) . At 
short times the plateau decreases with decreasing amplitude. 
By way of contrast, at long times the MSD is nonmonotonic 
with A and has a maximum at A* ss 1. On intermediate time- 
scales superdiffusive particle motion develops when A ~ A* . 

expect these additional dissipative forces to reduce parti- 
cle motion; friction immobilizes particles by making them 
stick to the surface. As we have seen, however, this is 
only true for the short-time dynamics, but not on longer 
time-scales, where particles even move super-diffusively. 

Given the polydispersity in masses, the Coulomb 
threshold is different for different particles. There is, 
therefore, a regime Aj ' \i s g £ [m s ,mi] w [0.52,1.43], 
where the driving force is just strong enough to mobilize 
the small particles but not yet strong enough to overcome 
the Coulomb threshold of the large particles. It is in this 
regime, that we observe the anomalous diffusion. 

At low enough densities, i.e. without inter-particle in- 
teractions, large particles would remain immobilized on 
the surface (Fig. [3] left ) . At the high densities that we 
study, large particles nevertheless move as they are kicked 
around by the mobile small particles. As a result the 
MSD is quite similar for large and small particles, as 
Fig. [3] shows. This is clearly a consequence of the lack 
of space at the high densities under consideration. No 
structural relaxation is possible when the large particles 
are kept frozen in place. 

The kicks temporarily lift the large particles above 
their Coulomb threshold. They can then undergo some 
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FIG. 2: The effective diffusivity D is nonmonotonous with 
a minimum in the glassy regime and a maximum at A* ~ 1. 
The temperature T is monotonously decreasing with decreas- 
ing amplitude. 
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FIG. 3: Oszillation amplitude of a single particle as a function 
of driving amplitude A (left). The threshold for mobilization 
depends on particle mass. At A — A* the large particles 
are immobilized on the surface. At the high volume fractions 
studied large particles nevertheless have to move, as they re- 
ceive constant kicks from the small particles. As a result their 
MSD (right) is quite similar to the one of the small particles. 

slip displacement before sticking again on the surface. 
One mechanism to generate supcrdiffusive motion is pos- 
itive temporal correlations in the kicks. To evidence 
these correlations we have to change the point of view 
and ask what the small, mobile particles are doing over 
time. Consider, therefore, the following situation, where 
a small test particle is confined to a cage-like surround- 
ing, which is taken to consist of only large particles. If 
these large particles are immobilized and not allowed to 
move, the test particle as driven by the external force, 
will move around and explore the available free volume 
(Fig.©. 

If we switch off the frictional interactions of the test 
particle with the underlying surface, the motion is quickly 
randomized by the collisions with the cage. As a result, 
a pore-pressure is built-up, which can act as a restoring 
force when the particle configuration starts to change. 
Consistently, in the full simulation with frictional inter- 
action between small particles and surface switched off, 
we observe particle localization and no supcrdiffusion. 




FIG. 4: A small particle surrounded by a fixed cage and driven 
by the periodic force explores its free volume. Without fric- 
tional interactions with the surface (left), the particle motion 
is randomized by the interactions with the cage walls. With 
friction (right), the particle can very efficiently synchronize 
with the driving and a periodic orbit results. This is possible 
because the particle can stick to the surface, when the driving 
force pushes the particle into the walls of the cage. There the 
particle remains until the force switches sign and drives the 
particle in the other direction. 



By way of contrast, with friction turned on, the test 
particle very efficiently synchronizes with the external 
force and explores only a small part of the available vol- 
ume. Instead of being sent off into a different direction 
during collisions, the particle sticks to the surface for as 
long as the force needs to switch sign and push it back. 
This synchronization is apparent in the phase of the oscil- 
lations of the particle. Without friction the particle coor- 
dinate is out-of phase with the forcing, just like a driven 
and undamped oscillator. With friction the particle ve- 
locity is in phase with the forcing, like an overdamped 
oscillator. 

In this situation a motion of the cage walls only weakly 
perturbs the trajectory of the test particle and, in partic- 
ular, does not lead to restoring forces that would oppose 
this motion. This is what allows small changes in the 
particle configuration during one cycle to be propagated 
and amplified in consecutive cycles. Thus, the kicks of 
the test particle on the surrounding cage are correlated 
on the time-scale of the structural integrity of the local 
free-volume. On this intermediate time-scale we there- 
fore expect super-diffusive, or even ballistic motion. On 
longer times the local surrounding has changed and any 
memory is lost. This is the onset of diffusion. It would be 
interesting to see, if on longer time-scales this diffusive 
process can lead to similar segregation processes as are 
expected for higher driving amplitudes [y, 0] ■ 

Inter-particle friction - Let us now investigate the case 
of finite inter-particle friction = 0.1). Particle-surface 
friction will now be switched off (p s = 0). The above in- 
terpretation of the origin of anomalous diffusion suggests 
that the role of friction is twofold. First, it acts as a dis- 
sipation channel that removes kinetic energy and makes 
particles stick. This alone would lead to localization and 
the arrest of particles. Second, friction acts as a non- 
linear perturbation and therefore allows for irreversible 
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ments, whenever interactions with an external container 
are to be avoided. 

In the experiments of Ref. 
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FIG. 5: MSD with interparticle friction and strong damp- 
ing f = 50 chosen such that velocities are damped on time- 
scales shorter than the driving period. The driving force 
(A — 50) is only applied to small particles. If all particles were 
driven, no rearrangements would occur (overdamped limit). 
Inset: example of a periodic particle trajectory when friction 
is switched off completely. The amplitude of the motion cor- 
responds to roughly 5% of the particle diameter. 



dynamics. These two aspects may be separated into dif- 
ferent processes, for example, by overdamping the mo- 
tion of the particles with a strong (linear) damping force 
-Fdamp = — C^i that acts as primary dissipation mecha- 
nism. Particle motion is now again in phase with the 
forcing. In this scenario, and without inter-particle fric- 
tion, the MSD (taken after full force cycles) is identical 
zero and particle trajectories are strictly periodic (as ex- 
emplified in the inset of Fig. [5} . The remaining role of 
friction is now to induce small irreversible perturbations 
on these periodic trajectories. Accordingly, and in line 
with our reasoning, this leads to super-diffusive parti- 
cle motion on intermediate time- and length-scales (see 
Fig E]). 

Conclusion - We have shown that, quite unintuitivly, 
frictional interactions can speed up particle motion and 
lead to superdiffusive dynamics. This contrasts with 
what friction is expected to do: slow down particle mo- 
tion by making them stick. We have rationalized this 
nontrivial behavior by considering friction as a weak irre- 
versible perturbation to genuinely reversible trajectories. 
This fluidization competes with the localization due to 
randomized particle motion. As a result we encounter a 
reentrance melting transition from a granular glass (lo- 
calized) to a fluid state upon lowering the temperature. 

These results may be important for experiments in sev- 
eral ways. First, note that the role of frictional interac- 
tions with the surface is quite subtle. Fluidization only 
happens on relatively long time-scales, while on short 
times, particle motion is completely normal. There, fric- 
tion only leads to a suppression of particle motion be- 
cause of the Coulomb threshold. Such effects would make 
it hard to evaluate the relevance of friction in experi- 



cle motion was indeed observed, however only at volume 
fractions around the critical close packing limit. In con- 
trast, we find anomalous dynamics for a range of densities 
(Fig. [5]). The high stiffness of the brass particles in the 
experiment makes a key difference with the simulation. 
We speculate that only at cj> c this high stiffness provides 
a sufficiently tight packing, such that frictional effects 
can compete with the randomization due to the driving. 
The role of friction would then again be to provide the 
nonlinear ingredient that can "rectify" the motion. Our 
analysis suggests that important additional insight can 
be obtained by changing the driving amplitude. First re- 
sults [2(J indeed suggest that the anomalous dynamics is 
enhanced for lower driving amplitudes. 

Finally, the transition between reversible and irre- 
versible particle motion is a nice example of a transition 
into an absor bing state, as described in the experiments 
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Interestingly, in our case the particles 
continue to interact in the absorbing state. This can, for 
example, be seen from the periodic trajectories of indi- 
vidual particles, which are complex loops and not just 
straight lines (inset of Figl5j). It remains to be seen if a 
volume fraction can be identified, which plays the role of 
the critical point of this non-equilibrium transition. We 
leave this question for future work. 

We acknowledge support by the Deutsche Forschungs- 
gemcinschaft, Emmy Noether program: He 6322/1-1. 
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